NPS ARCHIVE 
1961 

BORTHWICK, R. 









'f® 

SOME DYNA-MIc' AIRFRAME STABILITY ' ' Vi' 
CONSIDERATIONS IN THE eONTROi:’ PROBLEM i’, ‘ 

ROBERT BRUCE BORTHWICK 



it :i 



. ’ I ' , 

ill, ''ll' ' 



,u n ' 1 , -i 




•: ■''■•'I ■ vfnt'k- 



IllililliiiB 







LIBRARY 

U.S. NAVAL. POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 



r^- 

.. 



■ 














i 




w 

4 






S'." 











t 




SOME DYNAMIC AIRFRAME STABILITY CONSIDERATIONS 



IN THE CONTROL PROBLEM 



ic ic ic if ic if 



Robert Bruce Borthwick 



SOME DYNAMIC AIRFRAME STABILITY CONSIDERATIONS 



IN THE CONTROL PROBLEM 



by 

Robert Bruce Borthwick 
Lieutenant, United States Navy 



Submitted in partial fulfillment of 
the requirements for the degree of 

MASTER OF SCIENCE 
IN 

ELECTRICAL ENGINEERING 



United States Naval Postgraduate School 
Monterey, California 



19 6 1 



^<oi 



V. S. NnvnI Post'^raduntc ScliooT 
Moutcroy, Caliiornirt 



SOME DYNAMIC AIRFRAME STABILITY CONSIDERATIONS 
IN THE CONTROL PROBLEM 
by 

Robert Bruce Borthwick 

This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 

ELECTRICAL ENGINEERING 
from the 

United States Naval Postgraduate School 



ABSTRACT 



The dynamic stability analysis of the airframe for consideration 
in the control problem was made considering six degrees of freedom. 
The intention was to explore a method and determine the extent that 
such a scheme could be used in a system known to be non-linear at 
larger disturbances . It was determined that a very workable system 
exists whereby differential equations may be simply converted into 
block diagrams and vice versa. The system is workable for various 
flight conditions. It can also handle disturbances up to one tenth 
radian in most cases but is restricted to linear approximations . 

The system is compatible to any multi-loop control problem 
and is not restricted to airframe control. A standardized form now 
exists upon which multi-loop compensation theory research may be 
conducted using either basic hardware or differential equations as 
the basic system. 

The graduate work, for which this thesis is a partial require- 
ment, was performed while the author was attending the U. S. Naval 
Postgraduate School during assignment to the Aeronautical Engineering 
curriculum (Avionics) . 
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1 . 



Introduction 



The types of motion that result from a disturbance in some 
equilibrium flight condition and the transient motion of the aircraft 
in response to control movements are analyzed in the study of the 
dynamic characteristics of an airplane. Over the course of years the 
effects of dynamic stability in aircraft have been known to designers. 
The mathematical computations were discussed by early aeronautics 
pioneers which included Lanchester, Bryant and Glauert. In spite of the 
fact that dynamic stability effects were known, most people paid little 
concern to this field. As design groups solved other problems of air- 
craft design, the dynamic stability effects were automatically solved as 
a secondary result. However, since World War II, aircraft development 
has pushed into the field of very high performance aircraft. In many 
cases the dynamic stability problem doesn’t solve itself as happened 
previously. Furthermore, the fire control problem has entered the scene 
which has generated a control problem of no small magnitude. 

This work will consider the aircraft in six degrees of freedom and 
explore paths of compensation and solution. It is the prime intention to 
explore the work of Chu^^^ in multi loop servo systems and determine the 
value of it as applied to the aircraft stability. 

The first section will set down a general definition of terms and 
the basic equations of motion as related to aircraft. It will further set 
down the basic assumptions upon which the equations are formed. 
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The equations motion for the airframe which are completely 



derived in Appendix I are: 

- - (y 4. y „ HV)p = y . nt- F-in 0 (z^l' 

V (V ^ R'J - P-On s Y + n:- cos ©sin O 

^<7. “ = T ni‘~ "OS O COG (3’='”^} 

L = '■■p 4 + (q „ P)"R p(-'’ » r2) „ + R) 

+ ’^PR ^ c K35 == b K^R 

¥. = + bK, + (u - C)P: + E(P^ ~ R'’) = F'-R + DPQ 

+ cKj^R - CK 3 P (I» 5 ) 

R s HR + c + (B - + F(C" - P^) + E(f:R „ f) 

- DPR. + bK2P - aK^: (1-6) 

■^e " f’e " = ^elescs = ^ex 

^ ""r * " >"rVcy = 

? Hg * ^ ^2 (R''- + P) (I»9) 



Fhc pny’lfr rnrl linear velocity of a free body In space arc 



P (>10) 

' O CCS 0 + ^ pin 0 cos ^ (’"^11' 

R = 0 cos sin ^ (■''”>11) 

d^’^» s 7 r-os © cos 4^ -f- V (sin 0 sin 0 cos - cos 0 sin ^ ) 
dt 

4 - PJ (cos 0 sin ©cos * sin ^ sin if* ) (T-11^ 

dyj :- = cos © S’’‘r' •+* + V (sin 0 sin ©sin ^ + cos ^coa^) 
dfc‘ 

■*■ (c 's 0 sin © sin tp - sin 0 o''s ^ ) ( J=>lii) 

^ _-T ^ sin 0 cos 0 + W COS0 cos © (T-1!)) 

dt 



2 



I 



I 







•* #■ 



* mm 






« 






* 



« 



** • 



I 




Consideration of the aerodynamic forces must be resolved into 



the effects they produce on various parts of the body in some flight 
condition and also the resultant moment effects they produce. 

The first or necessary condition is to establish a reference 
condition and proceed from there. A reference condition is some flight 
condition where equilibrium exists in all equations and the summation 
of forces is zero. This is a steady state condition. 

The set of axis upon which this rigid body is oriented is commonly 
called the Eulerian Axis system. General aerodynamics usage modifies 
the z cixis so the positive direction points downward . 



Fifi’.re i“l 
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Crifdn ?t 
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The first rotation moves about the Z axis in a direction ■'f and is known 
as yaw. The next is the movement about the Y axis which is term © 
and known as pitch. The third motion about an axis is roll about the X 



axis and is maneuvered in terms of ^ . Caution must be used in de- 
termining these positions because the order in using this system is not 
commutative . 

The derivation for the yaw rate R , the roll rate P and the pitch 
rate Q are shown in detail in Appendix Two . They are: 



P = ^ - (j) sin © 


(1-16) 


Q = ©cos 4^ sin ^ cos & 


(1-17) 


R = 4^cos ^ cos© - © sin ^ 
and where small angles are involved 


(1-18) 


P = 0 


(I- 19) 


Q = © 


(1-20) 


R ^ ^ 


(1-21) 


Similarly, Etkin shows: 

© = Q cos - R sin ^ 


(1-22) 


^ = P + Q sin ^tan © + R cos ^ tan © 


(1-23) 


= (Q sin^ + R cos ^ ) sec © 


(1-24) 



The general equations of motion have to be applied to the aircraft 
motion . The effects of axis movement upon the weight force of the 
aircraft will be shown as well as the reference state condition. It 
must be remembered that the gravity force is always oriented toward 
the center of the earth. The reference body axis may or may not be 
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oriented in this direction. 



From the discussion of the Eulerian angle concept, it is 
determined that a reasonable reference condition would be 



F = O = 

X 


X - mg sin © 
o o 




(1-25) 


Fy = Y - 


- W cos © si 1 0 


o 


(I-2G) 


o 


o ^ 




F^ =0 = 


Zo cos cos 


5^0 


(T-27) 



This allows for a steady state condition where the aircraft is not 
accelerating but neither is it necessarily in straight and level flight. 
From the previous illustration of Eulerian angles coupled with the deri- 
vations in Appendix II, the gravity forces can be determined as shown 
below. 



CCS G cos 


CCS Osin ^ 


-sin^ 




ny,. 


cos 4^ ©ran ^ 


cos'll cos ^ 


cos &sin^ 


xng = 


Wy 


-san^cos <p 


+ sin'r^ sin^sjr 0 








cos4^ sinCl^cos 0 


sin'f’sin ©ccs 0 


cos©cos^^ 




'ny^ 


- sin4^ s'c (p 


-cos sin ^ 




! 





Discussion and derivation of this general cosine matrix are 
found in several sources (2,3). The overall result of the equations oi 
motion is shown below. 
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G *0^^^ Or-r0 

^ J --r ' r= Y ^ ; ( r i r> O ^ V s n O f; \n 0 p ^ n r> or ^ ) 

'^Y (cos r-'^ n cos y^coc ^ © *=: ^ 

r ' (cos cos^ (cos © rjri ^ ) ( < 

(’^^^I + F7 <* ^"0 ^ (r5n©^) (cos ^cjn Geos 0 « s-'r^s'jo^) 

+ p"" (c"'<5©^ S'' (s3 o n ©c'^s^ « oos^s"n^) 

+ (cos©^ cos^q) ( cos ©cos^) 






L - ?K^ P ~ + c - b^^P 

M ^ ^ * FQ ('-0) + ’^. (P^ => R-) + -Kt^ ^ cK, P 

j. 3 



(:~=‘u) 



f r s 



- np - - K^ 



r>'^ ( Vi^ 



) + ^ (^R ~ P'1 - hKp^ - ~ K-l' 



C^-33 



The above equations are complete except for the external forces 
X, Y, and Z and the external moments L, M, and N. These forces and 
moments are composed of aerodynamic and propulsion effects and the 
moments due to control surfaces. These equations are very general in 
nature and contain the following assumptions for simplification from the 
most general and complex case. These assumptions are; 

1. The earth is assumed to be fixed in space and the Earth’s 

atmosphere is assumed to be fixed with respect to the earth. 

2. The airframe is assumed to be a rigid body. 
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3. The mass of the body is assumed constant during the time 



of the problem. 

4. For an aircraft or missile, the XZ plane is assumed to have 
mirror symmetry. 

The above equations then cover the six degrees of freedom found in an 
aircraft or missile. Coupled with the other auxiliary equations covering 
roll, pitch and yaw effects along with control motions, solution is pro- 
bably possible. To date, this has been accomplished either by a machine 
using specific coefficients for various terms or several simplifications „ 
Discussion of this will be conducted in succeeding sections. 
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2 . 



Consideration of External Forces and Moments on a Body in 



Motion 

In the previous section the equations of motion were discussed. 

The aerodynamic thrust and control effects were simply grouped in X, 

Y, Z, L, M and N. It is the aim of this chapter to discuss these in 

detail with respect to the effect on the dynamic control stability problem. 

Consider the term X when the aircraft is in a disturbed condition. 

The reference condition is X and thus X = X + AX for some disturbed 

o o 

condition. The question is then to determine the A X components. 

The force and moment variables can be expressed in a Taylor's 
Series^ . 




For practical reasons the effects of second and higher order terms are 
omitted from further considerations as will be discussed later. Con- 
sidering the motion of the aircraft itself the X forces would be 




‘ff-r ‘(W *(«)’ 
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Similarly 




For convenience all terms will be expressed in letter subscript form„ 

It can readily be seen then that considering X, Y, Z, L, M and N and 

the effects of U, V, W, P, Q, and R plus the acceleration terms U , 

/■ 

V, W, P, Q, and R there are 72 derivative terms required. This does 
not include the control effects which can add another 36 derivatives. 
Needless to say, considerable effort has been expended in justifying 
simplifying assumptions. In addition to this obtaining reliable co- 
efficients for various derivative terms is not always possible . 

It is not the intention of this work to cover the computational 
methods for these various derivatives. The scope is too large, A 
search of all authoritative books p3, 4,5, ^on the subject yields a 
table of terms which can be reasonably determined. This appears as 
Table I. These considerations are mentioned to bring out the awareness 
of the different causes and effects of external forces . There are con- 
siderable gaps in the table . Determination of these to complete the 
table would be a topic of a complete study in itself. Much work remains 
to be done in this field. 

Furthermore, some of the methods of computing the known 
stability derivatives are questionable. Some computations from theory 
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prove quite accurate. Others are very poor and empirical methods or 
model testing are used. Even these fail at times and '"educated gue --. e, 
from previous airplanes are used as a basis of determining stability 
derivatives for a new design. Tail effects are the most notable in thi 
last category. 

In all aerodynamics problems , considerable effort is made to 
non-dimensionalize the equations. The advantages of doing this are 
of considerable merit. For the most part, it is rather simply done. 

It has been found in both theory and actual practice that 



F =4/®V SC 
2. F 



II-4 



Where S is a reference area, usually the wing area, V is the absolute 
reference velocity and Cp is some non dimensional coefficient. C^. 
varies with Mach Number and Reynolds Number depending on the force 
considered, (lift, drag, etc.) is the density of the fluid through 
which the body is passing. Thus a term such as Zw has the following 
relationship: 




Here again V and S are the velocity and wing area terms respectively, 
The same concept applies to the moment terms except that a moment 
arm term must also appear in the equation. A generalized example is: 






m 



ii-e 



In this example M is some moment and C 
coefficient to satisfy the equation. , 



, is a needed dlmensionles . 
M 

V and S are the same as 
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mentioned above whiled is the moment arm. 

Likewise in a stability moment the derivative can e 

form. 
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Another external force to be considered is that of thrust. 
This becomes an important factor in most flight conditions. The 
location of the thrust axis is seldom coincident with any reference 
axis in the body. 







m -r-*, •* • 
1 



'»'! TT’’ 



'V- 
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Assuming the thrust line lies in the xz plane, then the thru h: 



forces can be resolved in this generalized manner, 
X = T cos 0 



II-8 



Z = T sin©„ 
T T 



11-3 



For a steady state condition the equations become 



^Ot = To cos e T 



II-IO 



11-11 



= To sin e T 



11-12 



Oip 



11-13 



while in a disturbed condition where the body axis remains fixed in 



the mass , 



= Tj oos e J 



11-14 



= Tj^ sin 0 



II- 1- 



Mt - 



and = To + AT 



II-IG 
11-17 

Assuming the air density remains constant so the engine thrust will 
not be affected 

i -44 

11-18 



bs 



RE' 



or 



. 1 c i 
d RE' ^ 
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These are the considerations for the external forces on a bod/ 
in motion. While the dimensionless coefficient has been used exten- 
sively by the cieronautical engineering field, it will not be used further 
in this work. The primary mission of this work is to investigate the 
control response and paths for possible compensation. 

The generalized equations of motion then appear; 



o 



; . A ’> -• p 11 t-T t.T ' « 






ic? 0 



S0H.U + T( 



p 11 rp‘* rpm 



’ =• k 



•vj f \T j r 



) + jn/ s'n ©,^(cos ‘f^sin &sin - 'inH^cos^) - n ■ • cf' • 0 in 
^ 7 4. ri..,7 + Y« 4 V. j. 4- t + yf V »■ V IT-' 

» -jJ ■ '}i 

c& 

c ' ' - ’ ) ^ ■ •''in6'^(coc l^ps?^ ^rcs ^ ^ pf ■ 

(. S Ir 0 ^) (-'in Sf^sin ©CCS ^ - cos^^sin “r ricocG^-- .0 

>- T _.n&m 









q’’ n 

'■ r. ^ 



11 * z « • 






'T. 






c n . h K^Pt “ L Y • L 

^ 0 V ^ 









. Y(P" 



< a 



c K P ® M’ u 



f n 



^ '•" V * - M» - ' *'»^o + p» Yi 

'r / o ' j. - 

~ P <^r:r z 

d ro- 



? • >1 4 'r „ i 'T' •, 

■^ 1 - '- 0.4 



^ „ iC 4 4 e('R=F) + bK P = r. Ki = P ^ ’ • N" 

^ J- O 

4- 
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Thus, it has been shown that the equations of motion have the 



above external forces which will affect the rigid mass in motion,. The 
assumptions made are: 

1. The earth is assum^ed to be fixed in space and the Earth's 
atmosphere is assumed to be fixed with respect to the 
Earth , 

2. The mass is a rigid body. 

3. The mass of the body is assumed constant during the time 
of the problem. 

4. For an aircraft or missile, the xz plane is assumed to have 
mirror symmetry. 

These assumptions were made in the first section, in addition the 
following assumptions and simplifications are made in considering the 
aerodynamic forces. 

1. The flow is assumed to be quasi-steady. 

2 . The effects of auxiliary aerodynamic equipment such as 
flaps, speed brakes, slats and spoilers are neglected. 

' 3. The thrust can be affected only by change of speed amd/or 

revolutions per minute. No change is considered for change 
in fuel mixture, manifold pressure, side slip or propeller 
fin effect. 

These assumptions were made for several reasons. In the first and 
third cases, absence of adequate information on many of the coefficients 
makes ignoring or setting them equal to zero mandatory. In all probability 
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these coefficients are very small. Assumption two covers special 
cases. While these items may have considerable effect when used, 
they are eliminated from the general case at this time because an air- 
craft probably uses them less than two percent of its flying timeo For 
the special cases, these effects may be quickly added. 

The work in these first two sections is developed primarily for 
background purposes. There is nothing new or startling in this. For 
further examination of this subject, the reader is referred to the works 
in Ref. 7 and 8 which closely parallel the development in this work. 

The aerodynamic force and moment terms carry a primed notation 
for ease in publication. In later sections these terms will be divided 
by mass or moment of inertia. These terms after division wall not 
carry the prime (“ ) notation . There are far more of the latter terms 
used. Thus, it is for ease of publication only and does not affect 
theory in any way . 
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3. 



The Airtrame Equations in a Standard Form 



This section will take the general equations developed in the 
previous chapters and simplify them. The simplification will extena 
to an aircraft in straight and level flight and subject to small pertar' 
bations . This is the standard method of dealing with the problem as 
found in present day texts. However, these equations will be used 
in block diagram form and in determinental arrays .as developed in 
Chu's work. The results will be compared with present day 
classical methods to determine validity of this determinantal concept 

The first simplification is the concept of the lateral system 
being independent of the longitudinal system. Therefore, the six 
equations of motion are divided into the two groups. The longi- 
tudinal group includes the directional motion in the X and Z directions 
and rotary motion about Y axis. The lateral group consists of direc- 
tional motion along the Y axis and rotary motion about the X and Z 
axis , The axis system referred to here is the body fixis system of 
the airplane . 

Secondly, the aircraft is initially in a steady state condition 
This is the reference state as mentioned previously. When the dis- 
turbance terms are all zero Equation 11-19 reduces to: 
r. •< Qr. o O’ PnVr ) m'’ sin ^ ^ = Xo -<■ cos © III-l 

However, it was stated that this situation was restricted to longi- 
tudinal or lateral systems. Thus for motion along the X body axis, 
a yaw motion Ro and a side motion Vo would be zero. Further it was 
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stated th^it the aircraft was in a steady state reference condition 
Therefore ^ would be zero. Thus equation Ill-i reduces to 
n (OoKo) + IT." s:' n © ^ ® X*c ^ To ccc & T III 

Similarly, it may be observed that equation II-20 for the steady 
state reduces to 

Yo = 0 III-3 

Likewise equation 11-21 reduces to 

■m cVo) TC cos © + To sin© T ® 0 III-4 

c o 

The moment equations all reduce to zero provided the rotor effects 
are neglected or considered zero. 

Lo = Mo - No - 0 III-5 

Recalling that only small disturbances were to be considered 
and these disturbances w'ere to be in only the symmetrical or longi- 
tudinal system, for the longitudinal group the derivative due to asy- 
mmetric variables of V, P, R, ^ and ^ are zero. In a similar manner 
when the asymmetric system effects are considered the symmetric dis- 
turbances are ignored. This leaves the variables V, W, Q, and 
out of the asymmetric equations. Finally, because of the small dis- 
turbances, the change angle ©will be small so the approximation 
cos 0^1 and sin 0 S' j0 is valid. 

The result of the previous discussion is subtracting equation.. 
III-2 from 11-19, eliminating the asymmetric changes and supplying 
the small disturbance considerations . The result is; 

PI (u ^ >4.) + mr © cos © „ * X + x„*.: + 7 * 'Y 

III- € 

+ Ty-' ccs©,pU f T^^ppf,.noEex'<S RH' 
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For the z direction taking III-4 from 11-21 and using the same 



assumptions' 








.,e- '• 

* ’ Tv^' S"' ©rp “ ,1 



ri r 

e 



For the moment in the y axis , 

^ 3 * ^ 3 

^ ^ Ki{l ^ McQ 



the following applies: 












-1 



III 8 



" T y,-H- 

To complete the case for the symmetrical group, the action of the 
elevators must be included. From appendix I the elevator hinge 
moment effects are reduced to: 

tg ^ ri"*^ ^ '“Tl. ^ ^ewW + Pe'' ' ■*■ 1 11 “9 

In most cases where moment of inertia terms are handled it 
is usually expedient to use an axis system which is the principal 
inertia axis. This means that the product of inertia terms are zero. 
For aerodynamic motion problems , this applies for some cases . 

This case is not one of them. What is gained in simplification of the 
rotary motion equations by use of principal axis is lost in attempting 
to handle the force equations with the wind velocity considerations 
in sine and cosine terms . By orienting the X axis into the reference 
wind velocity the result is: 

O O 0 

C 111-10 



where U = + u . 

Therefore, in the steady state reference condition \/ = U and 

CO 

Vq = = 0. With the assumption that u, v and w were very small, 

their products and squares can be neglected. 
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and since U ^ is much greater than 2 then V ^ U and tie, i, 

o ^ C O 

terns are interchangeable in this section. 

As previously stated , the reference situation is a steady stare 

motion and includes first order effects only. The orientation of the 

axis U is a finite value while = W =0. Furthermore, the author 
o O o 

is hard pressed to visualize an aircraft or missile tumbling along with 

a Qp, motion although P exists in spinning missiles and bullet i . 

o 

Further X “ X = Z * = X>7 “ X-vi = T _ = 0 as an assumption, 

w q ^ ‘ ^ S 

These terms are usually very small. 

The "standard" determinant form as stated in Chapter 3 of 
Chu^ is: 



1 - P"ll 


FTp 


FT 

r. PI 


•^3 


^31 




1 ■*- PFp 


^22 


R^3 


PF 30 


• P^13 






1 ^ PF 3 


p- 

33 



where the FF indicates a performance function. The system may huve 
/y-i^ nodes. This example contains only three nodes. The "standard" 
block diagram for this form is: 
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175 3-1 STATIDARr* ^LOCK DT''.AR/..M DPAA'A'T FROT 3x3 DETFRMBIAOT 



21 



The nodes are labeled in Roman Numerals,. The basic rules 



of nodal signals allow several inputs, but only one output. This 
output from the node goes directly to a performance function with a 
single subscript number. Signals may be picked off from the node 
output and fed forward to another node . An example of this is container 
in the signal going through PF,,,. In this case the signal is also multi- 
plied by the function PFjg and then fed to node III. Feedback takes 
place in a similar manner. The signal output from the first system FF^ 

Jim 

is picked off and fed through PFjj^ back to node I. Similarly the other 
nodes function in the same manner. 

For feedback blocks , the first subscript number denotes the out 
put signal source for the feed back and the second subscript nurober 
the input node to which the feedback signal is going. The notable 
points about the determinant are: 

1. The main diagonal always has terms of one plus the 
direct path performance function times the self-loop 
feedback performance function . 



H- PFjPF^j) 



2 . All terms above and to the right of the main diagonal 
are feedback terms . 

3. All terms below and to the left of the main diagonal oi 
the "standard" determinant form are of a feed forward 
nature. All of these feed forward terms are negative. 

The term immediately below the main diagonal term is the 
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oirect path performance function « Thus, where the sec or, d 
column main diagonal spot has 1 + PF^PF^^, the direct pat! 
term PF^ w*ill appear in the location immediately below .r. 
-PF... 

4 . 



These rules are strictly mechanical. For a rigorous derivation of this 
system the reader is refen'ed directly to Chu’s work. 

In Chapter 9 of Etkin , there exists a two degree of freedom 
situation where ©and W vary and the speed in the forward direction Is 
constant. Tnis is the short period approximation. The equations usir g 
the LaPlace transforms and nondimen ional coefficients are; 









cC 



c< 






a 









- n 



Ill-i. 






n*^r 



) & ’n 



III-12 






hri 





III-? 3 



This author also grouped his acceleration terms such that G*t 

'-i 

included C> .. and . Similar situations exist for 

and G . Note that ^ and ^ ~ for this consideratlor' . Fro^n-i 

these equations a block diagram is evolved. 
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Fn-vrr T^LOGK 07 TJO D^^-iREUS O'- FREf:;?;: Ri THE LC' T^n'T/'L 

SYS?"' SHOE?’ BY ElV ’nr' 



Further, Etkin derives: 



"6,, ■ i . C,. „ •S-.t 



III-14 



^ ^ , C , ^ „ fj ^ ^ *1 >2 C; ~ I** r v J « ^ r /V ? r» 

i -yi^ ^Os' ■■ 



V.. 



'=<>1 






n '(>= 






m-i; 






Cttu ) 



rSTT *"''’-7 ')« 



Not it is intended to use Chu's work involving the use of 
"standard" block diagrams and "standard" determinant form to prove 
that the 'standard" system applies for the airframe as well as a 
multiloop ser/o system. 
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The Z, M and equations of motion with the applied as,. urr<p- 



tions and S =0 are: 
o 



t /> ■ ^ T ' 



,7 hn ■ « 



10r^\ - » 10_;j . 10 



1C 



a Fe 






in-17 



III -18 



Recalling that co = Uo t-'n ^ o0.and for small angles 
a,' A and 0 q , 

The equations in LaPlace notation then become: 

=. (slCtJ * (slO^; * lO^P (slC':)©-AFe III-13 

(s'- = (sM-rJ + + Mo() tjoAK ^ C III- 20 



vS 



UoAA 



(sU^ 





III- 



These equations in standard determinant form now become: 
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A < 


© 


1 - ('slG-^ ^ 10 -n ) 


-i" slOot +• 13s< ) 
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1 = 


“ s^TTq 4 7o) 
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« ( ^ M ) 
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Ordinarily, it is not necessary to divide by s or Hov, 



ever, it is done here to demonstrate more closely the standard lou^ 

To take the differential equations and arrange them by rowo n*.. 
columns to see the actual physical system is to no avail. The nr inge 
ment does not exist in the above form. It will be shown how to achie’i.e 
this form later. 

Consider the "standard form" array as set forth above. The 
block diagram for this in Chu's "standard form" is: 
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The above block diagram then is a three node system but the outputs 
from the first and second main path blocks are nothing that is le iJily 
recognized in an airframe. Because of the fact that the main block 
following the third node is 1 , the output here is 0 . 

A few simple block diagram manipulations clears the situation 
and shows true signal outputs that are physically realized in the di- 
frame. Block diagram manipulations allow the changing of a pick oft 
point provided the transfer function is corrected. Further, two block, 
can be combined into one equivalent transfer function. 





Adding a dummy node this now changes to: 




'"f' ' rp 
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The result of csin.j these two maneuvers rearranges the blocV u: ' 

to the following ^ 




i’in'.re 3“'‘5 
Or ^IC-IT'IO 



R”'ARP.'‘"''-tT’! RLOCK C7 FirrTTf.E 3-3, 

■5.J 2 --^ 



USP'Ti ’.A:'' 'ULi 



Maneuvering the above diagram location , the block diagram takes the 
form of Fig. 3-7, which coincides with Fig. 3-2. 
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T' '0 ■ :"p- '!^ I ''I '-Y'lT'*' 

^ ‘‘ of '*ic-’^re ^-6 to cor.prtre ■’^'^e 3-*^^ 
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It is noted that this block diagram approaches the form set 



down in Etkin except for nondimensionalizing the terms. The dasheu 
area in the diagram is what Etkin labels and G o( 

.-The determinant from the above dashed portion of the block diagram 
is: 
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Evaluating for the transfer functions and , the results 

will be in the form 

Cl} f 1 ^1} 

This form is the type expressed by Chu and Thaler ^ , The A 

term, which is the whole determinant, reduces to: 
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+ Moi,) 
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=s(U^ * ■’r) j 
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sI'L } 



■(Uo 



Zq KsM0<, 

,9^ 




s(Uo 



7 N(, 
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MoO 



III-2. 
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To consioer the^o^ term the cofactor A. y oi this deterndnant is 



obtained by crossing through the first row and the second column. 




The main feed forward transfer function between the input signal fror 
node 2 to the output is +1 in this case. This is also true of node 4 
input to output. Hence = +1. 

The transfer function 

^ - ____ . 










^ + s(Un ^ m:-i 



I II -24 



(A 



)(s' 



sI'U ) 



.(To * 



It can be noted that this equation III-24 Is the same equation as 
III-15 given by Etkin save the non-dimensionalizing terms. 
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To ev.Tlaale the transfer function & i 7] , the same niethod is u, ed 
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Therefore: 




1'^ Al'i 

A 




.s -'"iK:-" 



. -o<,) ■>- (gi-^o C •" TXVgZ i ^ Z-^ ) 
sK ) ^ rJ'r ^ V + ^’c<) 
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lU^ 



Here again, excepting the non-dimensionalizing , the result is the 
same. Ec^uation III-26 is equal to equation III~14 as given by Etkin 
The system of "standard" forms then works for the airfrarrre 
differential equations after some manipulation. It can be noted that 
there is a method of making up this standard form without the necessity 
of the manipulation which will yield physically realizable output signals 
When the block diagram of Fig. 3-6 is converted into a determinant 
for solution of the transfer function, there appears a series of alterna-' 
ting rows where the equation “1+1 = 0 appears » It also appears that 
these simple equations occur in a sequence where the plus one term is 
placed in the 2n! row and column, n is the number of differential 
equations in the problem. The minus one term is placed in the ?n!! row., 
and (2nl“l) columns. All feedforward terms (those below' and to tne 
left of the main diagonal) are moved one column to the right and reriain 
in the same row. This doubles the size of the determinant. However, 
there are several zero terms introduced and the actual work invol' xd in 
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evaluacion of t>'e determinant is not markedly increased. 

The next item is to look at the lateral motion problem, P Is 
slightly more complex because there are five differential equations 
of motion insteaa of four „ 

Equation II-2 with simplifications becomes: 
fT , T’or) ^ pirn- “ YrR ■« J = 0 III-27 



' £ • r) / 



rr ■'OS “ Y z*?" Y p „ Y 
0 p 






0 



V 



o 0 



The moment equation II'='22 becomes 

Jj ~ X-v^r — *0 



I> TT p « T, ^ L ^ U 






and equation II-24 becomes 



R - TP » ^ NpP » N^P . ^ ^ N^' J 
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The aileron hinge motion equations from Appendix I are 



^ ^ ll^P - ll^r - 11^ ^ “ 11^ 3 - A 
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And also from Appendix I, the rudder hinge motion is 



- - ir 



,P = I’r-^ 



^ - l->v ^ - A 



m-r;. 



34 




t 



Fron these terms a determinantal array can be made. It sns i 

be noted that this time a set of equations of -1 + 1 = 0 will be inoertei.* 
in the system. This permits direct drawing of the block diagran and 
omitting the need of manipulation. Table 3-1 then presents this on .y 
From this Fig. 3-8 is drawn directly. 
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TABLE 3-1 DETERMINANT ARRAY OF LATERAL SYSTEM EQUATIONS FOR SMALL D3STRUBANCE: 
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The resulting block diagram very clearly parrdlelr Tha’ o: ti.e 
lateral sva>'‘.em shown by Etkin in Chapter The '' -n:antl'ird"'' foriT, tc 

achieve this block diagram differs. 

There is, then, the question of why does this differ from tue 
original standard" form set forth by Chu. 

Chapter Two of Chu^^^ discusses at some length the idea oi 
arrangement of varying performance functions into a "standard"’ 6100* 
diagrara. The individual pieces of hardware are at hand with known 
individual performance functions. The orderly arrangement oi these 
permit the construction of a "standard" determinant for the whole system 
From this , a multiloop servo system can be handled in generalised 
results. In the case of the airframe, the determinant formed by the 
differential equations was known and the block diagrara had to be 
formed. Admittedly a block diagram (Fig. 3-3) was made up but the 
output signals are not the physically realizable type that are gener ally 
known to pilots and engineers. The same performance function results 
from the system. Both are mathematically correct. 

The difference lies in the fact that Chu is taking individu il 
pieces of hardware where the functions are known. The blocks tre 
made up and then the "standard" determinant form is made. In the 
of the airplane, the equations of motion are known and a determinant i - 
immediately formed. Each column contains the performance functions 
multiplied by the direct function signal output term. Ch'u's standurd 
form dictate^, that all terms above the main diagonal are multiplied tv 
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the node output or Input to the direct function. Steps rust be taien 
then to change this in the differential equations. The addition of ire 
alternating -I -r i = 0 equations solves this very nicely and in nc v. ly 
alters the value of the equation, 

A series of evaluations of determinant forms is given m Appen iix 
III, The results show that the addition of these terms in no v/ay chkin-ie^ 
the final value of a determinant or a cofactor value , 

Fefening to Fig, 3-6, there are other block forms that may chan ,e 
the values within the blocks, but will not alter the final outcome. Con 
sider the first node in Fig. 3-6 as an example. The four forms expresse 
in Fig, 3-3 are equivalent. The original form is shown in Fia, 3-3(a). 

The self-loop feedback path is combined with the direct path in Fig, 
3-9{b). The s^ term is eliminated and the performance functions are then 
in the form shown in Fig, 3-9(c). A further change can be made by 
dividing the nodal input terms by the denominator of the direct path 
feed-forward function. The result is then in the form showm in fig, 
3-9(d}. The performance functions in each block differ somew'hai, but 
the final result remains the same. These manipulations allow the user 
more freedom in the forms of the performance functions. It allows 
changes in form to fit the needs of the particular problem. 

The pick-off point for the feed forward terms does differ in 
requirements. It is seen in Figs. 3-9(a) and (d) that the pick-off can 
be either off the output signal of node I or II, The final result will not 
change. However, this is not the case in Figs, 3-9(b) and (c) , The 
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piC!' *01: point must be off the output signal of node II The ;in il 



result will difier if the output of node I was tapped In the Irtter c tse 





3“9(c) 



Figure 3-9^’h) 
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For con. erting a set of linenr differential equations into tne 



standard block diagram form, the following rules cm be applied, 

1, Between each differential equation, an equation -j. + 1 = 
■shall be inserted. An additional equation -1+1=0 snaii 
be added at the bottom of the set so the number of ;.ddi 
tional equations shall equal the number of differentinl 
equations originally in the problem. 

2. The arrangement of the values of these added equations in 
the determinant shall be such that the positive one value 
is always in the main diagonal block. The minus one term 
shall always be placed in the block immediately to the left 
of the main diagonal term. 

3o All other blocks in the rows of these added equations shall 
be of zero value. 

4„ All finite performance functions of the origin, il different! d 
equations which are located below the main diagonal shall 
be moved one place to the right. The vacated blocks ohah 
be replaced with a zero value. This last rule is optional 
to the user. It does save the possibility of error whict’ 
could appear because of cases shown in Fig. 3-’3(b) and 
(d) . 

With the simplified rules of making a standard block diagr irr« 
from the standard determinant, a valuable tool is developed for nviny 
purposes. It c m provide a student with a means to draw" rand vfsu »ll;e 
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a CO” plex feedback system such as an airframe from a group of difieien 
tial equ^itions. For the research groups it makes mrilo j cor?.putei 
progr \fr,nlng a simplified matter with respect to .implifier usage . In 
this latter case, pitfalls do exist. The airframe is no simple mattei' 

The reader is referred to Howe'^®^ for an excellent coverage of the dr 
fraire problem for analog computers. Most important, though, is the 
fact that this standard form of Chu”s is valid for any m,alti-loop linear 
feedback system whether it is electrical , mechanical or a combination 
of the two. This then puts a regular form to the multi-loop control 
systerr, from which design or synthesis can proceed in an orderly rriinner 
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4„ The Control Problem for an Airframe 

The previous three sections have discussed only the nifi iix.e 
itself. It is quite reasonable to say that while these point out m my 
of the responses, they achieve little toward solving the control probJe" 
In simplified terms, the control problem is resolved into the block aia- 
gram of Fig. 4- 1 , 




--'Tgur;.' 



SHTLI’^'C: OCAITT'OL SYST’""" 



Indeed this diagram is an oversimplification of the aircraft 
control problem. However, it does point up the fact that the previous 
three chapters have dealt only with the controlled system, block. Any 
complete study of the airframe control problem, must Include the other 
block components of the amplifier and the error detector. The probler 
is a nultiloop type and of no small magnitude. 

For the control problem, the diagramming becomes a simple 
extension of the block diagram developed in Section III. Consider the 
standard determinant for the longitudinal system. 
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The deTernlnantal array for the longitudinal system is is follow o 
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From this a block diagram can be made up from the deter min«.ni d 



array . 




OP AI?'’RA*-1F ’-TTH TKF~^ PiiGPK'f^ 



This diagram is of the three degrees of freedom longitudinal 
case. Combining this to a control system, the following diagram 
results . 
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Figure U-3 CaiPLETED LONGITUDINAL AIRFRAlffi SYSTEM RESPONSE 



Thus the control problem for the simple attitude & is 
presentevi in Fig. 4-2 and 3. The two represent the different point of 
view of the aeronautical and the electrical engineer although the uy ter 
is the same. In Fig. 4-2 the aeronautical engineer is faced with tnc 
primary problem of compensation within the airfrairie. For the electrical 
engineer who inherits the fixed airframe , the attention would shill to 
Fig. 4-3. In this figure, Gj would represent the pilot response in movin 
the control lever. G 2 would represent the servo system response of u 
hydraulic or electric system. Gg would represent the error detection 
and correction signal gained from an electronic element of some Find 
Similar situations of block diagramming would exist for the change^ in 
speed and altitude. For the most part, the compensation schemes to 
achieve stability or improve stability are fixed with respect to pati.s. 

The chanvge in aerodynamic coefficients will vary the paths in the air- 
frame response group. If the electronic devices are to be used, then 
the pick offs must be at the system outputs and the signals fed bacf 
to near or at the input or reference signal node point. 

In Fig. 4-3 the pilot constitutes the outermost loop. Hov/evur, 
this is not mandatory. In a fully automatic system, the pilot circuit 
is opened . Few paths exist where both aerodynamic or electr onic 
means could be used to achieve stable control. One can be noted on 
Fig. 4-2, That is elevator position response 7|^ which could be altered 
aerodynanically by changing parts of the present feedback block or 
adding another electronic control block. It is then fed into the fore e 
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input noae marked /X Fe. There is no guarantee that this p;-.th 
any stability problem albeit the flutter problem may be iided,, It moielv 
serves as an example of two possible types of compensation over the 
same path. 

1'he lateral systems are not shown or discussed in this chapter 
but the same principles of control systems apply. It should be rertem 
bered that the airframe response shown here covers the small pertur- 
bation response concept and thus is greatly simplified. The following 
sections will deal with expansion of the concept for some workable 
systems . 
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5. 



Lxpansion of the Stxind^ird Form for Greater Distuibanres 



The system developed by Chu was for linear function^, It. ij 
the purpose of this section to explore the equations of motion and 
determine whether the equations of motion can be used for lai ^ei ois 
turbances than considered in Section III. If so, then to w’hat degree 
can this systerr be used. 

For all practical purposes, the small angle approximation hold^ 
quite w'ell up to one-tenth radian . Therefore , S ^ sin 0 and 
cos0:g: i. With this simplification in mind, Equation II-Iu is then 
rewritten . 

n fu + 'W RV) + mp sinO^ = (cosSq sin^^) + 



e Ceos©,, cos ) * Xn + X*U + X’W x«V X«e + X“ ^ 

o mo' o U ^J- Q 



X'lj X-V + T cos© + T cos A rnU + Tr L ms 
I V o T RFK * RPK 



V - ■ 



The thrust forces are combined into the forward velocity motion 
u and a steady state reference condition of 
n - RJ^) ^ mg sin©Q - X^ + cos © ^ V-L 



is assumed. 



Subtracting V-2 from V-1, the following equation results 

m 'u » •* W - ow =■ R V r « vr) 

cw cq ov o 

“ V’ m'- (cos©^ sin^^) + ©-mg (cosO^ cos_^^) - X.,' 




4 




X^».V = 0 



/ 3 
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The product tern qw is dropped as a small value procuct situation 



The sane applies to the vr term. 

Converting to La Place notation and dividing by the v-'-iSi 
followinQ equation will result 

^ PCS © cosine » X 1' ' (sX', * X-,)A < = pXc Q 

o / o u ^ \f. 

- - X >^)7J ^ J^/3 « 0 

The y and z direction force equations will form up in the same m 

" v' o« ^ 0^0 ^'=< - «o SZ5 “ 

- g (C0,i6^ 

V' 

- - pYj. (y^)4 - (sY ^ + Y^)^ =0 

•* g coE0^, s-’'r^Q ^ “ g cosa^ f'osp 6 ^ “ '^^,U 

= (s<' * 7,,) Ac<. ». sZ^^e - + 7 ^)r( = 0 



tt.e 



4 
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The moment equations can be formed in much the same manner From 



Section Two, Equation 11-22 would convert to this equation: 

- . (sl,| ♦ L4 'I ' 



In a similar manner Equations 11-23 and 11-24 can be brought to a 
similar maximum linear condition. They are shown in the following 
equations: 

+ sQq^) + 3T - sPo</^) 

V-8 



f *• ^ Iff sPq0 - 



Also the control hinge moment equations 111-9, 111-31 and 
111-32 are recalled in their same form. 
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These nine equations represent the linear equations of motion 



of the airframe and include all cross-coupling effects. The_e equatlor. . 
still retain the requirement of a linear approximation. The equation ^ 
with the c.uthor's modification of adding the equations -1+1 = 0 ir 'i e 
up Table V-I, This is an 18 x 18 determinant. From this table the 
block diagram of Fig, 5-1 can be drawn 
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TABLE 5-1 



DETERMINANT ARRAY OF AIRFRAME WITH LINEAR APPROXIMATIONS FOR SIX DEGREES 



OF FREHXM 






I 



I 
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u his biocK diagram then represents the complete lirir ime 



response m iintaintng linear approximations . Even Then th> Is ysTe'" 
h 1 .S liTatatlons . As an example, consider some airplane llying uitr. 
a 3te-;.':y state wing angle of attack of 15 degrees, (assiine a ‘'con 
vention" wing and the reference axis is the zero lift line) Any pilot 
will quickly testify that an increase of 5 degrees angle of attack at 
that point will probably result in a rather wildly oscillating maneuver 
popularly known as a spin. This condition is generally outside of the 
range of the aerodynamic forces and moments usually computed . Sor^e 
aerodynamic coefficients such as lift (0^) are decidedly non-linear in 
this region. Therefore, the restrictions imposed on this derK'ation of 
terms must fnclude a restriction that the airframe is not operating si; 
the very edge of its flight envelope. 

The system has been diagrammed directly from Table ¥”1„ Tne 
terms in Column One are all connected at node I of Fig. 5-1. Likewl^se, 
the terms in Column Two of Table V-i are all connected to not*e 11. dxtr i 
nodes appear in the beginning and end of the block diagraRi. These ire 
inserted as provisions for control functions. These control nodeo are 
lettered alphabetically. The block diagram contains the required loop 
for aileron and elevator control. The rudder control loop can be e ^ily 
inserted w’hen desired. 

Exarr.inution of Fig. 5-1 reveals many interesting points here- 
tofore known but sometimes only by intuitive reasoning. The concept 
of separate systems for the symmetric axis and asyrr.metric &xi > 
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S'yoie'* s is . er/ well shown. The level flight steady st ite condlticn 



usea m Section Shiee stdteJ that 

Po = Qo = Ro = Vo = Wo = 0 V 

Sxarimng the paths between nodes VIII and IX re’ceal th it ill 
of the perfoirtance functions are zero under these specified condition^ 
of steady state flight. There is no direct, feedforward or feerib^sCi 
paths between nodes VIII and IX. Thus the restricted case of steady 
state level flight conditions with striall perturbation is clearly divide! 
into ti e two systems The longitudinal system includes the forwar ; 
velocity, vertical velocity, pitching motion and elevator position. 

The bloch diagram for the system can be further restricted to two 
degrees of freedom by holding the X velocity change, u, at zero. 

This aliov/3 ch tnge in vertical motion and pitch. This would then be 
Figs , 3-6 and 3“7 of Section III . 

The so-called lateral airframe response system comprises tie 
nodes from. IX to X'/III, The rudder and aileron motion, sideslip, y iw 
and roll comprise the motion for this sytem. With the steady state 
level flight conditions previously mentioned, the block diagram, is 
identical to Fig. 3-8 of Section III. 

An interesting note of these separate systems is the cros^ 
coupling termu that begin to »arise with larger disturbances. An ir 
portant one is the pitch due to sideslip. The effects are labeiev. 
in the diagram. Several aircraft are subject to this. Usually, tbij 
is y, nicVing or nose down motion that appears as side -lip angle ( /V 



;ncrt;i OS. In the Sirr.piifieu systems this efiect v..^3 nejiLectel It 
is ol ne tli-Mnle v.ilae lor very sm ill disturbances , but o in i.eco’ e 
notice ible by the time the sideslip angle ) re'aicher:. five .ejree 

7ne ellect exists and is shown by this expanded systeru, 

ISxceilent rrsethods of analysis are available lor these >t;and:mi 

forrts brought forth in this work. Methods of analyzing outputo in 

(l 

response to one or more inputs are available. Work by Thaler con 
tiins £ very useable system. 

S the output signal at the n th node 

A 

~ characteristic determinant 

^ the cofactor of the determinant wdth ",-i” 
the node location receiving the input 
and "n” the node location of the output 
3ign<al . 

■e the performance function on the input signal 
prior to entering the input node , 

A 

r. the input signal 

s the direct path performrince function betv.een 
the output node, and the output signal. 

Thus a single input signal would yield an output signal 




For two or more input signals additional subscripts would be used aid 
the functions are additive. The rules of superposition of linear ourputs 



justify this method. 




Where |V 
A 
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This e <.u it*on vi-ould give the resulting output of one or noic ‘nput^ 
at virio..s nodes in the system. 

The legrees of freedom can be restricted also. In the bloci. 
dii ir.U'-i f O'igo ^^-2) consider a cut made in the path at point s". 




i' lire 



OJil I 



I ’1 .^< 7 * 



-.y 






« -rn y* Tfp T 
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rt .3 opening in the path now eliminates the signJ to t'*e a ire it 



pith, t^e teedb ic)^ and the feedforward path. The result of tfis n 
determinism is the insertion of zeros in all the perfonriisnce funciions 
in t ,.e color n. It leaves only the constant, one, in the main d? »gon d 
bloch. The leader will recall that the main diagonal has the tern 
1 a Phh,PF,^^ 5 . « This is shown below as opening the output path from , 

C' ^ o 

third node. 
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The output path from node three is opened at point s and the deter= 
minant is changed in this manner, 

A special case occurs for the diagram illustrated in Fig ' , 
The second node has only one input signal. Because the path v,as cut 
at point "s", the output from the second node v/ill also be zero Ttius 
all feeciforward values from this node will be zero. The direct and 
feedback paths of the second node are non-existant to begin with md 
remain so. Tfe determinant form, would appear as now shewn 
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Simple manipulation now reduces the six by six determinant to a 
four by four determinant. Therefore, when it is desired to consider 
a problem with one or more degrees of freedom held at zero, the simple 
method of opening the output signal path from the node and placing a 
zero for all the performance functions in the column leaving only the 
constant value one in the main diagonal point accomplishes this 
Referring to Section II , the reader can observe the general 
equations of motion for mass particles and examine one such as 
Equation 11-22 which is repeated below. 









(r 



-B) = E (PQ + R) + c « bK.F 



- !.«,/ + LJiv L«p + L'l’ + 



» 

It 4 



V-13 



Note the third term QR is a produce of (Q + q) (R ’ r ). The pro- 

o o 

duct when multiplied is • Q^r R^ + qr. From the equations 



of Section I: 

Q - S CCS ^ ^ ."in 0 C06 0 V-j 4 
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" c o'c ^ r s & ^ © £?:in ^ 

When the sirnll angle approximation is no longer valid, the hrindhn.. 
of this product term becomes rather involved Furthermore, the aero 
dynamic coefficients of the Lp, etc terms lose their linear 
approximations over large changes in flight conditions. Most aero- 
dynar^-ic coefficients are of a decidedly non-linear nature over the 
full span of an aircraft's maneuvering capability. The coefficients 
are definitely non-linear with Mach number when the velocity r -.nge 
covers the subsonic, transonic and supersonic spectrum. The reader 
can observe the general equations of motion in Section II. Considering 
the complexities of incorporating equations V-14 and V-15 for large 
movements and the non linearities of the aerodynamic coefficients , 
it is not long before one realizes that the non linearities are numero;-o 
and cover all the equations . The prospect of block diagramming for 
the purpose of placing all non linearities in one block to be h'lndled by 
a describing function scheme is hopeless. 

The possibility of doing the problem with an analog computer 
and using no approximations is fair. Undoubtedly the task, if atte'^pted 
will be of considerable magnitude. 

Thus, it is seen that this system of standard determinant anays 
and standard block diagrams can be used to handle the airframe probler' 
as long as a linear approximation is maintained. The system is flexible 
and can carry varied conditions of multiple inputs. It can also be noted 
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thit for the six degrees of freedom problem over a large variation, 
it is of United value except for diagramming a computer problem . 
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Prospect and Retrospect 



The airfrarr.e rriotion problem has been derived quite ^'enei llv 
in the first two sections of this work. The attempt was made to miinUun 
as many effects as possible in the general equations of motion. Tfis 
was done so th at the basic equations of motion could be used in i 
variety of ?ilrframes varying from missiles to ground effect iri&chines. 

In the succeeding sections the author has confined his discussions to 
equations of basic airplane types. Rotor effects were disregarded and 
thrust effects were combined with drag to determine the effects along 
the X axis system. This should not imply that the equations of motion 
cannot be used for other airframe types. They can in all cases where 
the assumptions of a rigid body, constant mass, etc. are still valid. 

The author has several thoughts on areas of further exploration 
He will discuss them briefly and mention any work of possible interest 
in connection with this field. 

The primary aim for control engineering In this field is the 
development of an orderly process for multiloop compensation. It 
shown clearly in Fig, 5-1 that several paths of possible compensation 
exist. What path or paths are the most successful and/or easiest to 
achieve multiloop compensation? Some work has been done in thi^« 
field. Specifically, an electronic compensation of a missile using thi . 
determinant method was done by Anderson and Roane, They confined 
their efforts to the roll and yaw coupling problem for a specific vehicle. 
Their efforts achieved a satisfactory solution to meeting the nlssile 
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specific’it.ions , but no definite method of analysis as to the best or 

poorest pt.th to compensate a m.ultiloop system was *jener t*e ! . A 

possible approach would be an examination of aerodynamic coefficients 

5 

and their known effects. In chapter ten of Perkins and Hage tnere is 
a discussion of the fact that drag or coefficient of drag markedly affect,. 
the damping in the longitudinal phugoid motion of an airplane. The 
airframe short period is affected largely by the elevator hinge moment 
coefficients C and C . An examination of these factors ana 

others night lead to a pattern which could successfully solve the 
problem for any multiloop system. 

The block diagram development should provide a basis for 
analog simulation of the linear problem expressed in sections III and 
V. The analog cortjputer arrangement may be further expanded to inclu le 
some of the non-linear effects using function generators. This probler . 
allowing six degrees of freedom, would overtax the capacity of the 
analog equipment presently available at this institution. If an unaiog 
computer v/Jth a 100 amplifier capacity and 50 function generators was 
provided, the problem could be nicely handled. Most assuredly, com- 
binations can be made by grouping various functions which will cut the 
above nu.mbers by one half. However, various outputs will not be 
available for analysis and many signals which are not physically 
realizable in the actual system will appear as outputs. 

The high speed digital computers could be utilized jn a project 
of statistical evaluation of compensation schemes. The cop'puter's 
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Cap ..tjlity of producing h'^ntireus of solutions per minute woula -iiiovv 
ex-.r.in aion ot u.it;i to determine trends and patterns for viiiuu-- co- 
penrutor schemes, A successful compensation criteria inijfct de 
derive-i in lids manner. 

An aerodv'namics yroup may examine four items of interest in 
this field. The control surface hinge moment equations used in tnis 
work were extremely simplified. Exploration of these effects to 
attain more refinement woul i be justified. 

In section II , the author pointed out the absence of conpu 
tational methods for various ■'tability derivatives. Further ^ the 

assu'tption of a quasl-stea ly state condition is made in compitinj 

(l 2) 

the known derivatives, A study by Etkin cites an improved sysre". 
to attain better accuracy of the aerodynamic stability derivates , 

Further work in this area appears justified. 

There is also a need to determine the effects of aeroelaslicity 
on the problem. To this end, instrumentation and flight testing an 
airframe night prove useful. The capability of the Naval Air Facility 
to maintain and operate many types of military aircraft should pem'iit 
use of ”00-^6 airplane or aircraft for this work. An additional use would 
checking the effects of automatic stabilization equipment in use on t*^^ e 
airframe. The HSS -1 helicopter has an excellent control system for 
this type of research. However, the equations of motion used in this 
work will require some refinement to control the varying inertia problOT- 



caused by the tilting rotor disk. 



A final topic foi consideration is that of signal flow theory 
This idea, which has been advanced in recent years, may be com 
patible to the same '* standardized" system as presented here for the 
block diagram method . 
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Conclusion 



It has oeen conclusively shown that Chn's 'stand ird bloc* 
diagran and deterninant method can be applied to a linear syste^n 
other than & servo control system. With some modification, any set 
of linear differential equations can be block diagrammed into a stand ird 
form for general analysis, Avery workable system exists for analysis 
of ca nulti-loop control system using one or more input signals .. The 
airfram.e stability problem can be examined for the linear approximation 
using this very effective system. 
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BASIC IQJAxIONS OF MOTION 
T- e ! ‘isic c oncept X of r;:ib5 in notion cind the eqa'Uions cJ 
motion ctn be derived in tector not ition. This is done usin.5 mas.. 



particles in motion and determining their effects in two separate iixl,3 



system?. The basic reference system shown in Pig. Al-I is the . t irre 
axis. This is the non^rot’iting reference system and is also known is 



the Galilean Axis. The axi? system is fixed in the rigid mass but may 
rotate about the fixed space axis system. It is the relative system. 

A special set of notation is used here in this appendix which differs 
from the body of this work. 

The definition of the special terms in this appendix is as 



follows. 

4 

i 

j 



the vectorial unit along the X axis 
the vectorial unit along the Y axis 



A 




4 



J 






k 



A 



0 






'4* 

0 

R 



4 

4 



4 



the vectorial unit along the Z axis 

'IP 

the vectorial unit along the X «axis 
the vectorial unit along the Y' axis 
the vectorial unit along the Z axis 
the origin of the relative axis system 
the origin of the space axis system 

the vectorial distance from the space axis origin to the 
relative axis origin 



r 



the vectorial distance from the space axis origin to 1 
mass particle in space 



!»iau..(lT!liiiii 




dm 



/I 



L 



A 

A 

A 



the vectorial distance from the relative axis origin to 
the mass particle 

the mass particle 

the vectorial distance from the mass center to the origin 
of the relative axis system 

rotational velocity vector of the relative axis system aboiit 
the space axis; an absolute velocity vector 

absolute rotational velocity vector of a rotor whose ixis is 
attached to XYZ axis system 



H 

G 

F 



A 

A 

4 

A 



/ 



momentum vector for the entire mass body 

moment of momentum vector with respect to the space axis 

moment vector with respect to the relative axis 

the summation of the external forces on a body 
A 

= the summation of the internal forces on a body 
This appendix wall derive the basic equations of motion for my 
mass moving in space. After the basic equations have been developed, 
the effects of control actions will be added. The final equations will 
then be converted to standard aeronautical terminology which is used in 
the body of the work for aircraft stability considerations. 

These equations have no doubt been derived several times in the 
past one hundred years. However, the author was unable to find any 
clear derivations in any of the reference material covered. Therefore 
the derivations in this Appendix are submitted for the reader's interest 
and for background information . 
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Consider 



r jX + ;n/' + lc7 



U’ 



JiL = i u! 



^ ■*■ ■•■ k'-o- 



Momentum vector for the whole 




dm ^ 
dt 



dr = fVo */ of rotor) 

dt 

L =J •■’ “ J ■*■ (u)x/^'' + rotor) 

hi ' fd-v. d'T" s r 

dt J J 



but summation of internal forces equals zero. 



Hence 



dl, 



f dv ^_r ~ T7 

dt J dv 
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For an airframe with a large rotor (ie helicopter, STOL aircraft or GEM> 
consider addition of the rotor to the airframe using the mass center oi 



the aircraft as the origin of the axis system. 






fi y (-H- X + /^ + 2/lx /^, 



9 I 

(it" 



c - 



7 



tr’^nsport rnlntive corjolis 

pccclrrntion ef^C'';t 
t'^ *9^ 



< ■ t *• 






C' r 



r . o 



-=5^ -:S> -=>■ 

A ” ^ (1 X ( n X /^ ) ^ /^+ 2 CL r./^ 



.n 



0-37 ^ ^ ^ ^ 

*• j-"-Vy + Xv-V? ' iX-(n.pVz -CL^Jy) + 

k^Kil^^'rr '■ /I2’'”) + 0 + 0 + 0'+0 vrl^re - 0 



(Cl, 



'A. '3 

















= 




k-:0- -- k-::-(V7 



-n. 



3'-'- 



^A^Yy 



’•A-’V)r O 
A;^V3)in ( 

“Ao'-'x)n 
^ J 



naplecting rotor effoct. 



Th'^re’^ore in scoD.rr for\T' 



?x = (V>: ''T -A^Vy)'!! 

■"y * ( / +A,'y -il^y.^)in 

dz ^ (Vz Y' -A^Vx)!?! 



Concurs with Ktk"‘ n 
h ,L( s3 

/•Iso Duncan P-^p Sh 
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is fixed in the rigid body 
ang\ilar velocity of but the mass 



, The rotor spins at an absolute 
center of the rotor system is at 



point Po Therefore j rc and/^c “ It is also assumed that 
_n. = constantc Therefore^ the additional terra^xCrix/Jg) 
enters the problem » 



UX X /"i. 



A X (3 jc C; = (A*C)3 - (A“B)C is a basic identity ^ 



i 


j 


Ic 


^1 


Mo 


K/3 


rc 

y: 


A. 





therefore 
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Note; coordinates i , j , k are those in the body and not the inertial 
or earth system. The values of-A.are primed for the sane 
reason. It is quite possible that will be 0 considering 



alignment of the rotor on the xz plane. 

F.„ - F, + F., - ir.,a + m.^a F, ® airframe F„ = rotor 

I 1 ^ i * 1 L 

In resolving this, it must be noted that the mass centers differ 

from those of the combined masses. 

It should be noted that the equations for the acceleration of a 

particle also contain rotational and relative velocity terms, however, 



Moment terms may be considered using Moment of Momentum. ■ 
From this the differential may be taken and the moment is then 



derived. Consider H to be the moment of momentum of the system with 
respect to the inertial (non“rotatlng axis) . The derivation will include 
a rotor spinning on the rigid body. 



and = 0 as these particles are in a rigid body. Further, the 
summation of particles about the mass center makes 




3 t 
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The angular motion of a rotor in motion , but attached to the body will be 






( U,\ C‘ 



-h ^'1 



- A 3 ) 



Henceforth this relative term will be noted as 



The equation for H is: 

P - /r dir. r ^ R + 

y RT 



K' (KAPPA) . 



Recalling that the momentum equation was 



L = 
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then 




The term undei the integral will be known as H. 

Recalling that - 7_ ■» /L • t . ' . ^ 

/f "T 1 

Yi=Jl^ - <■ t > j % 

This can further be resolved into the moment of momentum about its 
own Oxyz axis and the velocity moment with respect to the inertial 
axis . 

- (M + m)/^ c X Vo *J ^:iLx/^ dn * f x (l( x/^) 

V* * ^ (JXx/^'I dr 




Then substituting these in. 







JL ■* - 





However , we also noted that 
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Therefore; 

4 4 4 4 

s 'c L *" r. ”■ r’l, 

c't d^ dt '•■'t 



Eqfaating the two results of rrHV- it is seen that 
+ G ^ ^ X L + 

^ dt (hv, dt 

Recall that . ■ . 

—S' — ^ 

^ Vo 
dt 



and that 

I. ® (M * m)(Vo t/i.2c/^ */^a 

vr.t,T ^ s d_ (M * !>’)/^ Vo * Ho F 
dt dt L 

- ^ -4 ^ 

d/^, r ’-o V 






dH * (in ^ I’) 
dt 



dt 



ot 



^ 4^ 

Ho ^ Ho , 
rel 



cowiHnf 
(\ 



-4 --4 ^ ^ 

Vo X (Vo ^Ax/^^) M * Vo V (Vo +Ax/ „ rotor) m 

C C ' c 



■" ‘ ^ y Vc * ^ A x/^^) X Vo)m + 

1—-^ "4^ 

- (?' <- m) X dVo + d _ Ho + hOj.g2~] 
dt dt L J 



- X Vo (/lx/" ) X Vo) I' 

An 



A further assumption is stated at this point. This is that the 

point of rotation of the rotor is fixed on the body and remains so ,, The 

shaft does not have nutation about some axis but remains fixed on a 

mount of the rigid body. Therefore, the of the rotating body is 

O/ 

zero. Also, logically, in a rigid body the mass does not shift hence 




is also zero. 
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In crossing a vector with itself Vo x Vo the result is zero. 



o X .0 



n 



n 



• 9 



*^^0 r '.'o ^ 0 



= j (^'pV3 - i 2 (V1V3 = 

* k (Vj_Vo » V 373 ) = 0 



Furthermore, crossing two dissimilar vectors in opposite order yields 
opposite signs and equal magnitude. 










i j k 

V 7 V 

123 

A A A 



7 (V2/>3 - *j(/^iV3 -ffj) 






X 



(A 



7’3 



j ?■ 






■ <^ 3 V ‘ 



Thus Vo X 



^ X 0 

/ C / c 



Incorporating these two facts, the equation for G can be simplified to 



the following: 

q = |Ho + 'AelJ ^ 



-j> ^ 

riFo ® ^iLx ho :'n a 

?nr 



"enaral e^n.-^tiono 
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Thus it is necessary to solve dFo to find G. 

dt 



Considering Ho = J / (u! x 
absolute velocity. 

Recalling an identity of vector algebra: 
' X xF) s =» T' »?)(; then 

A - (A 3 )/ 



which accounts for the 



Ho 



(x^ * ^ - 'r. . ^ 



r 



y -^2 + j u 5 j + * 7^2 ^ 

t + z“)u-S =. xyuip ~ x?u)^J 



(in iy ■*■ kr) 



+ ij^x^ 2 ^)u) 2 » - xyldlj 

+ k jTx- '*• y^}u)^ - xz u>^ - 



The above terms can be recognized in part as moment and product of 
inertia terms . 



Ix 




] + 7*/ 


)dn 


ly 


A 




''Ir’ 


Iz 


L 


j ix ^ y 


'd’-.l 


Ixy 


L 

y 


f Mex^ 




lyz 


A 


f fy’/dii 




Ixz 
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Therefore: 

Ho - i - Ixzu)^) ' J («3>yi^^ 4 Ty 

•+ k {i 7 - ‘■' I:-: 7 u^ - ly? i>2 ) 

Considering the effects of the aircraft as a rigid body first 

^ ^ t ^ "n 

d"'r> e Ho ■^Jlx Ho end for a ririd body u3 = -t I- 

" - * — ^ 

i (L'iILt * I:'y/l2 * Ixs./l^) J *- I7/L2 ^77^!-^) 



t, ** 



♦ k (^Jxyil^ * ^ l 7 .A^} ^ 



i j 

Ai yi2 ^ 

( Tv=Xx]'''*T':2 ) ( ^ 



For any ordinary aircraft or missile, the product of inertia terms Ixy 
and lyz can be considered zero because of the mirror symmetry aspect, 
Ixz may be zero if the principal axis is selected. Thus with no loss in 
generality, consider Ixy = Ixz = lyz = 0 

Ll? ® - TyA 

^ i? |j'’A 3 4 T,rA^.a2 “ 
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However, including the effects of Ixz, the product of inertia term, 



which is often not zero: 

« Txr,/l3 » JyzA^A^ 
jX‘A^ “ ^'2^1 + lyA-i/lj - 



TxzAA “ 

T^rA^/lo Ixz/l^n^^ 



Using standard aerodynamic terminology, the following items are 



defined: 




Ix 

ly 



c 






Iz 




P-Ai 

Q^fi, 

^ r, 

R - 



Thus 

G - trhere system is chosen sx'.ch thrt ?? o 

dt" 

? = 7 1^ yP- ^ - F,PO + CR (G-B)J 
f 7 + Ki (;-c) ♦ E (P~ = 

+ 7 Ipr - E? + PO (B-0) + 
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These equations concur with Duncan, Page 66, and Etkin, Chapter 4 
To consider the rotor effects it should be rernembered that w- 
Jl. f wheie K. is the relative motion of a rotor on the body. Thus 
for a rotor there will be a combination of terms considering the absolute 
angular velocity of the OlCfZ system plus the angular velocity of the 
rotor with respect to the body. 



It should be noted , that for a symmetrical rotor or disk the 
product of inertia terms are zero where the disk is at least a three 
bladed propeller. The two bladed propeller has a product of inertia 
dependent on its angular position. This case will not be considered, 
but the three or more bladed propeller situation will be used in all 
considerations . Hence with a dynamically and statically balanced 



The moment of inertia of the disk must be considered from the 
axis of the rigid body. Therefore, the moment of inertia must include 



Thus 




dt 



rotor = 0 and the product of inertia is zero. 



2 

the md consideration. 
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The Hioment of inertia terms for the rotors will be thus define ! in the 
most general sense 

a - Ix rotor - Ix + nid^ and aligned with the reference axis 

b - ly rotor = ly + and aligned with the reference axis 

A ^ 2 

c = Iz rotor = Iz + md^ and aligned with the reference axis 
G = dH dHo "* dH.^pi 

Recalling the derivation of moments for the rigid body without tlte 
rotor, the same method applies. 
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'it-, . ^/i.. p,, 



' "K ■' ■' (^ - ■< kc 1^ ;, ‘^(-' Tr,-( 



f ^T 



■i.rr ■- 



Combining this to a system consisting of a rigid body plus a rotor, 
the following results occur 



^ 4 -> r 



r 

-3 . i>-^a 2 w. 1 x 2 aJ , 

ixz.a^ zyR^n^ 4 ixA,-a ► ly^n 









*tl - 3^it, * Xc tf 



k 



i 2 
j 

jRi A2 A3 
al^ b{f2 



cf(^ 



Combining and using aerodynamic notation. 





[.f 






■R » clcoQ ~ 


'1 


[b 5 


- bK , . (; 


=C)PR + E(p 2 . 


r’) ^R » cl^ 


1 




%* 








[CR 


KP ♦ ^ 


(P»A)I^ -» D 


■R . blf,P . 



This result concurs with Etkin, page 116. 

Therefore, using a standard aerodynamic notation the moment: 
G = iL + jM-t-kN 

and the scalar quantities are as follows 
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Thus, the force and noment equations derived in this appendix 
are applicable to a rigid body of constant mass where the body is in 
motion. It also considers the effects of a three bladed or more rotor 
system v\hich is attached on some fixed point oi the rigid body. The 
rotor effecixt may be included or discarded depending on the magnitude 
of rotor inertia and the relative velocity of the rotor. 

It should be noted that the derivations are done in a general 
form., but in the concluding equations, the usual aerodynamic notation 
is substituted to bring the equations to this specific problem. Further-' 
more, missile or aircraft shape and density has been considered sym- 
metrical about the x -2 plane. Thus, the product of inertia terms Ixy and 
Iy 2 have been considered zero. No generality is lost by this and if, 
for some reason, the body is sufficiently un symmetrical , the terms c m 
be quicVly added in 

4 

The following is taken from Chapter 4 of Et in . The concept 
of motion on the control surfaces is shown in the follov^ing 







8G 




i 








The LaGr ‘in ;e equ-iiion of motion in a moving reierence Irame is 






r 




c? 







Si 



.J 



4 




kinetic energv oi the system relative to the chosen frame of 
reference 

= generalized force 




W = work dona on the system by the external forces which act upon 
it 



generalized coordinate 
H = aerodynamic moments 
S ~ rotation of control surface In radians 
Kinetic energy T ' t ' j __T - J-S 

3 

For elevator motion 

J • HeJ /o) - "eg'JV) 

h 

Fe - generalized elevator control force 
JWi - work done by inertia forces 
Work done by the control forces is 

* • d<^T 

^ J ^ 

e^r t.c co'^itrol. S‘ r'^^ce 

rr't-'o'" or ct" hydraulrc to 



h 




b7 



I 

I 



I. 




I 





Evaluating W1 

consider cFl = d r cim 

dt. 



recalling from previous derivation that . 

"“P' cT? 

:?2i ^ ^ </L^ P * '' P 

dt” dt 



Then 
dFxi = 






X 



'* =. C Ry = x(R^ + R^) ^ y (PC' =• R) + 2 ‘'FR-"')| ufj 



dt 



dFyi and dp 2 are similar in fom„ 

Looking at an elevator assembly, it can be seen that assuming a 
lamina in the xy-plane, the displacement in the direction is 



amottnt 



The work done by the inertia forces 



cf ( cT e)fk>i 



cfWi ®rd7nn/^edm 
jTjiT y dt ‘ 







edm 



R'" 




n 



The Coriolis terns 2 F>,- and -2 Qx are zero as there is no work done 
in those directions. Furtheinore z = 0 because of the assumed lamina 



PP 



cz // 

dtT“ 



e dm ^ 



-■ P <’FP. = 



F,. pi : 'n 't O' cf t’-''. p"! px'-atcr 

V 
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By sy-^ j’ietn 



ar: 



cT. “ “ ^ 

oT cT«) 

V- - '"c*Vc„ 

sji 



r-, . (F’-Q) 



^ (PR-'--) * + F^ 

f y ^ 0 0 



•'. " V‘e'3c, "■ ^ex ^PR-Q) ^ H ^ F„ (Ftkrln R.P, 13 ) 



Similarly ior the rudder 



-X 



/■’r'^r ” (R + P^) Prx ” (R'^ “ P) P^,^^ = + F, 



and for the aileron 



' • ..^. (pr * ?' - 2 p 
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APPENDIX II 



ROLL, YAW AND PITCH 

The derivation of the yaw, roll and pitch rates is given here to 
show the actual motion taking place. Frequently the motion is con- 
sidered small and small angle app'^oximations are used. These will be 
noted at the conclusion. 

First, consider the Eulerian Axis used in the aeronautical field 
shown below. The rotation takes place in the order of yaw ( ) , 

pitch ( ^ , and roll ( ^ ) . 




Ylgoxs AII-1 EUmilAN AlIIS SYSTEM SHO'^.THG RCTATION ABOUT THE 
Zj, Is AI® X AlilS 
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The initi il conditions it X , Y and Z and the relations of 

o o o 

f , - and ^ can be rol sted to P, Q and R of the body axis at 

X, y , z. j'n the above diagrir. consider the angular velocity F which 
is about the OX txis. It rnust be the sum of the velocities about the 



OZ^, OY., , and OX, axis respectively. 

Tf 



OZq maPes an angle ® with OX. The rate thus 



Y CCS f. 

& COc 



ll ^ Q.) si n O- OY is perpendicular with OX so 

^ I 



M 



~ 0 and OX is parallel with OX or the motion 



2 



Thus the summation P = ^ ^ sin ©• 



To evaluate Q the angle YOZ^ must be determined. Specifically the 

0> 

value Y cos YOZ is desired. 

cos " cos YCZ.^ cos 

(. C O 

- cos (TT - ) cos 0 

® sin p cos © 



Thus the contribution of OZ^ axis motion is sin ^ cos © . 

The angle between OY and OY^ is ^ so the motion here results as 
0 COS p . OY is perpendicular to OX^ hence the angular motion 
0 cos = 0 . 

K> 

Thus Q = ^ cos ^ ^ sin ^ cos © 

In the third angular velocity motion R the yaw motion appears as 
cos p cos 0 while the motion of pitch is © GOS 
The rolling term vanishes w cos _H1 = 0. 

Thus R = - Q oin ^ ^ cos ^ cos © 




4 



Hence the relationship of the body yaw, pitch and roll rates 
(p , Q, and R) to the designated reference system angular velocities 
, Q and 0 is; 

F = ^ 3in B 

a> 

Q = 0 cos ^ 4"^ sin ^ cos © 

R = ^ cos ^ cos © - © sin ^ 

These results concur with Duncan page 81 and Etkln page ilG 
It should be noted that for small angles the approximation belovv 
is often used; 

Q = © 

R = f 

A second consideration which must be made regarding movement 

about the reference axis is the velocity or displacement with respect to 

time. It is desired to obtain the velocity with respect to the fixed frir e 

or inertial axis x, j*y and 

Consider first that in scalar values 

dx* = U ^ 

TT” i 




dz^ = W 
dt 



It should be noted, however, that in the three position system shown 



on the previous page that 






t 





V.J ~ V cos ^ W in ^ 

W« = V' sin ^ + W cos ^ 
an.i W = “U^j sin O f W-> cos ^ 

U„, = cos A + sin O 

".y 

and U, =U,. cos f' - V „ sin ^ 

V| = U 2 sin - V., cos 

= w^, 

Coirblning the 1 , 2 and 3 position terms 

dx* = U C0.3 © cos *4^ + V (sin <p sin © cos ~ cos 0 sin ^ } 
dt 

1 - W (cos 0 sin © cos + sin 0 sin 0 ) 

dy^ “ U cos Q sin ^ + V (sin 0 sin © sin ^ + cos^ cos (f } 
cit 

+ W (cos 0 sin © sin ^ ~ sin 0 cos ^ ) 

= -U iin © + V sin^ cos © + W cos ^ cos © 
dt 

Thus, the velocity or rite of change with respect to the inertial 
or earth’s reference axis in this case is expressed in these terms as 
derived. These velocity terms are derived in Etkin Page 102 and coincide 
with the derivation shown. 

Therefore, this appendix accounts for the changes in anguLir and 
directional displacement of a rigid body with respect to a fixed reference 



J 3 



axis system. 
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APPENDIX m 



NOILS ON DEIiLkMINANT MANIPULATION 



The bocy of this projecr cont lins several determinant exp insions 
where the order is doubled. These manipulations do not alter the vdues 
of the deterninant and the resulting values of the performance functions 
are the s am e . 



The simple detennin nr 



Al 


hi 


Cl 


A2 




O; 


Aj 


B3 


C3 



evaluated will yield the follovving reoult 
A1B2C 3 + BiC2’A3 + CIA2B3 ' 
•CIB2A? - AiG2B3 - BIATCS 



A3-1 



This expansion by the equations -1 + 1 = 0 will result in the following 
determinant . 



Ai 0 B1 

I 0 

0 A2 B2 

0 0 -i 

0 A3 G 

0 0 3 



0 Cl 

0 0 

0 C2 

1 0 

B3 C3 

0 -1 



0 

0 



0 



0 



A 



5 



The resulting evaluation wall yield the same value as result (A3“l) 
Further- ore , if the terms below the main diagonal are allow'ed to 
remain in their original columns in lieu of shifting one place to the 



rivihi i .one in the n- ,in bociy the resulting A i>-> still the er e 
Ir.e results oi eualitritlng cofactors lor pfirticular perfomiince 
functions rerr -lin unch mjedo l or example it is desired to evaluate the 
cofactor A of the oricrinal deterntlnantc Thus; 
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■Ali — 


“ -hiF- ' 








A2 


1 

m 

1 


C2 


= (-1) /A2C3 “ 


- AZC:p A"* 


A3 


BIS 


C3 







In the expanded detenrilnant the equivalent cofactor is 



- 


M = 


- “ -o- 


1 


e- - 


■ “C±“ — " 


-H3- 










■s 


’o 

1 


0 


0 


0 








0 


A2 


m 


0 


C2 


0 




A 


J 

j 


1 




1 










iL V* 




0 


0 


1! 

-1 

1 


i 


0 


0 








0 


A3 


B13 


0 


C3 


0 








0 


0 


1 

P 

h 


0 


-1 


1 


1 






1 “ 1 


i 


II 

0 


0 1 




A2 


0 


C2l 


- 


'■J 


A^ 


0 


C2 


= 


0 


1 


0 




0 


0 


T 

A. 


0. 




A3 


0 


C3 




fl' 

1 0 


A3 


0 


C3 1 








II 

1 


/A2C3 - 


ascj7 










A' - 


The result of Eq^uatlon A3-3 equal 


s Eq'uation A3- 2 . 


With 


the value 



the main dragonal left in their original columns, the follov/ing occurs; 






6^ 



-0 — — C-^' — — O' 



io 



!0 



■no 



0 



A3 



0 



C3 



0 



0 0 0 0 -1 1 

I 

By Inspection of columns two and six, two reductions are made to the 
following; 

0 C2 I 
0 1 0 



A2 



A3 



0 



C3 



A _ = C-l) /A2C3 - A3C^/ A3-4 

1 S 

The resulting equation (A3-4) is the same as equation (A3-2) , 

A similar examination using a four by four determinant system 
was mace. This was done to insure that the simplified effects of a three 
by three or lower order system did not cover up an error in a four by four 
or higher order determinant. The same manipulations can be performed 
and the resultant value of the determinant is not altered. 

Thus, the individual has the option of placing his values in 
either column. The location of the values is dependent on desired 
location of the pick-off point on the block diagram. The final result 



is not altered . 



